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e'/e in the standard model 

M. Fabbrichcsi'' 

^INFN, Sczionc di Trieste and Scuola Intcrnazionale di Studi Supcriori Avanzati, 
via Beirut 4, 1-34014 Trieste, Italy. 

I discuss the estimate of the CP-violating ratio e'/e by stressing the role played by the chiral quark model in 
predicting the experiment and in showing that the same dynamical mechanism at work in the A/ = 1/2 rule also 
explains the larger value obtained for e'/e in this model with respect to other estimates. 



1. Introduction 



The CP-violating ratio e' je (for a review see, 
e.g., |l]-D) is computed as 



' 2 e ReAo 



Ho - - Hz 



.(1) 



where. At = ytdVts is a combination of Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements, 



Ho = ^^V2/,Rc(g.)n (l-f^^+^O , (2) 

cos On 

n2 - ^-Vy. Re(g.)2, (3) 

COS (52 ^-^ 

and (Qi)/ = (27r,/|Qi|X). The AS* = 1 effective 
lagrangian Cw is given by 



Lw = -^Y^ Ctifi) Qiifi) , 



where 



G 



dilj) = ~/^^^d. K; 2:j(/i) + T y,{n) 



(4) 



(5) 



In (|5|); Gi? is the Fermi coupling, the functions 
Zi{fi) and yi(/x) are the Wilson coefficients and 
T = -14dVt*/KdK*s- According to the stan- 
dard parameterization of the CKM matrix, in 
order to determine e'/e , we only need to con- 
sider the yi(/i) components, which control the 
CP-violating part of the lagrangian. The co- 
efficients yiifJ.), and Zi{fj,) contain all the de- 
pendence of short-distance physics, and depend 
on the t, W, b, c masses, the intrinsic QCD scale 
Aqcd, the 75-scheme used in the regularization 



and the renormalization scale fi. They are known 
to the next-to-leading order (NLO) order in as 
and ae [Q. 

The Qi in eq. (Q) are the effective four-quark 
operators obtained in the standard model by in- 
tegrating out the vector bosons and the heavy 
quarks t, b and c. A convenient and by now stan- 
dard basis includes the following ten operators: 

Ql = (SaU/3)v-A (^/3C?q)v-A ' 
isd)y^A T,q (w)vtA ' 



O2 
Q3,5 

Q7,9 



(6) 



i (sc;)v-A E, e, (gg) 



V±A 



2 {Sadi3)y_j^ J2q (^qiqpqa)v±A ■ 



where a, (3 denote color indices (a, /3 = 1, . . . , Nc) 
and iq are the quark charges (e„ = 2/3, id = 
Cs = —1/3). Color indices for the color singlet 
operators are omitted. The labels iy ± A) refer 
to the Dirac structure 7^(1 ± 75). 

Notice the explicit presence of the final-state- 
interaction (FSI) phases 5i in eqs. (||) and (||). 
Their presence is a consequence of writing the 
absolute values of the amplitudes in term of their 
dispersive parts. 

Finally, ^n+rj' is the isospin breaking (for 7tj„ 7^ 
md) contribution of the mixing of tt with 77 and 

^'- . 

Experimentally the ratio e'/e is extracted, by 

collecting K]^ and Ks decays into pairs of 7r° and 

TT ; from the relation 

2 ^1 

^^ ~ 1 + 6 Re - , (7) 

■qm e 



where 



?7oo 



and 



V+ 



{ttOtt"\£w\Ks) 



{7r+Tr-\Cw\K_ 



(8) 



(9) 



(tt+tt \Cw\Ks) 

The announcement last year of the prehminary 
result from the KTeV collaboration (FNAL) g 
based on data collected in 1996-97, and the 
present result from the NA48 collaboration 
(CERN), based on data collected in 1997 and 
1998 0, settle the long-standing issue of the pres- 
ence of direct CP violation in kaon decays. How- 
ever, a clearcut determination the actual value of 
e' at the precision of a few parts in 10^ must wait 
for further statistics and scrutiny of the experi- 
mental systematics. By computing the average 
among the two 1992 experiments (NA31 |0] and 
E731 H) and the KTeV and NA48 data we obtain 

Re £/£' = (1.9±0.46) X 10"\ (10) 

where the error has been inflated according to the 
Particle Data Group procedure (tj ^ tr x \/x^/3) , 
to be used when averaging over experimental 
data — in our case four sets — with substantially 
different central values. 

The value in eq. (lOh can be considered the cur- 
rent experimental result. Such a result will be 
tested within the next year by the full data anal- 
ysis from KTcV and NA48 and (hopefully) the 
first data from KLOE at DA$NE (Frascati); at 
that time, the experimental uncertainty will be 
reduced to a few parts in 10'*. 

2. Hadronic matrix elements 

If there is no sizable cancellation between the 
relevant effective operators, the order of magni- 
tude of e'/e is bound to be of the order of lO^''. 
A simple argument for this is presented in ||]. 
The problem is that any cancellation, or the lack 
thereof, among the operators heavily depends on 
the size of the hadronic matrix elements and there 
is no available estimate of them that is free of 
hard-to-control assumptions. To estimate the 
hadronic matrix elements in a systematic man- 
ner without having first to solve QCD (not even 



by lattice simulation) we needed a model that is 
simple enough to understand its dynamics and, 
at the same time, not too simple so as to still in- 
clude the relevant physics. We chose the chiral 
quark model (xQM) |l^ in which all coefficients 
of the relevant chiral lagrangian are parameter- 
ized in terms of just three parameters: the quark 
and gluon condensates, and the quark constituent 
mass. The model makes possible a complete es- 
timate of all matrix elements, it includes non- 
factorizable effects, chiral corrections and final- 
state interaction, all of which we thought to be 
relevant. 

The xQM |ll| is an effective quark model of 
QCD which can be derived in the framework of 
the extended Nambu- Jona-Lasinio model of chiral 
symmetry breaking (for a review, see, e.g., [p2[). 
In the xQM 3'ii effective interaction between the 
M, d, s quarks and the meson octet is introduced 
via the term 



C 



xQM 



-Af(9«I](ZL + giSt 



qR 



(11) 



which is added to an effective low-energy QCD 
lagrangian whose dynamical degrees of freedom 
are the u, d, s quarks propagating in a soft gluon 
background. The matrix E in (0) is the same 
as that used in chiral perturbation theory and 
it contains the pseudo-scalar meson multiplct. 
The quantity M is interpreted as the constituent 
quark mass in mesons (current quark masses are 
also included in the effective lagrangian). 

The hadronic matrix elements are matched 
with the NLO Wilson coefficients at the scale 
A^ ~ 0.8 (~ nip) as the best compromise be- 
tween the range of validity of chiral perturbation 
and that of strong coupling expansion. 

3. The fit of the A/ = 1/2 rule 

In order to assign the values of the model- 
dependent parameters M, (qq) and {asGG/ir) , 
we consider the CP-conserving amplitudes in the 
A/ = 1/2 selection rule oi K —^ tttt decays. In 
practice, we compute the amplitudes 
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A2 = -^KdK;-^E^'(/^)I^°(Q'(^))2 



V2' 



-LO Af) ^ri+iY 



(12) 



within the xQM approach and vary the parame- 
ters in order to reproduce their experimental val- 
ues. 

This procedure combines a model for low- 
energy QCD — which allow us to compute all 
hadronic matrix elements in terms of a few ba- 
sic parameters — with the phenomenological de- 
termination of such parameters. In this way, some 
shortcomings of such a naive model (in particular, 
the matching between long- and short-distance 
components) are absorbed in the phenomenolog- 
ical fit. As a check, we eventually verify the sta- 
bility of the computed observables against renor- 
malization scale and scheme dependence. 

The fit of the CP-conserving involves the de- 
termination of the FSI phases. The absorptive 
components of the hadronic matrix elements ap- 
pear when chiral loops are included. In our ap- 
proach the direct determination of the rescatter- 
ing phases gives at O(p^) 5q ~ 20" and 82 — —12". 
Although these results show features which are in 
qualitative agreement with the phases extracted 
from pion-nucleon scattering the deviation from 
the experimental data is sizeable, especially in the 
/ = component. On the other hand, at O(p^) 
the absorptive parts of the amplitudes arc deter- 
mined only at 0{p^) and disagreement with the 
measured phases should be expected. As a mat- 
ter of fact, the authors of ref. |I3] find that at 
0{p^) the absorptive part of the hadronic matrix 
elements are substantially modified to give val- 
ues of the rescattering phases quite close to those 
extracted from pion-nucleon scattering. At the 
same time the 0{p^) corrections to the disper- 
sive part of the hadronic matrix elements are very 
small. This result corroborates our ansatz |1^] of 
trusting the real parts of the 0{p*) matrix ele- 
ments while inputting the experimental values of 
the rescattering phases in all parts of our anal- 
ysis, which amounts to taking cos(5o ~ 0.8 and 
cos (52 ~ 1. 

Hadronic matrix elements in the xQM depend 
on the 75-scheme utilized ^^. Their dependence 
partially cancels that of the short-distance NLO 
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Figure 1. Dependence of ^o and A2 on {qq), 
(GG), A^^D and M at ^i = 0.8 GeV. The gray 
and black sets of lines correspond to the extreme 
values of Aqcd and M. The length of the lines 
represents the effect of varying (qq), while keeping 
all other parameters fixed. The vertical spread 
corresponds to varying {agGG/ir), with the cen- 
tral line corresponding to the central value of 
{GG). The gray area denotes the region spanned 
by varying all the parameters without correla- 
tions in a ± 30% box around the experimental 
values of Aq and A2 given by the cross hairs. 
The figure shows the HV and NDR results, cor- 
responding to varying (qq) (GG), Agrri and M 



in the ranges given in eqs. (|13[)-(14) 



(GC 



Wilson coefficients. Because this compensation is 
only numerical, and not analytical, we take it as 
part of our phenomenological approach. A formal 
75-scheme matching can only come from a model 
more complete than the xQ^- Nevertheless, the 
result, as shown in Fig. |^, is rather convincing. 
By taking 



A 



(4) 
QCD 



340 ± 40 MeV 



and fitting at the scale fi = 0.8 GeV the ampli- 
tudes in eqs. ([l2|)-([l2|) to their experimental val- 
ues, allowing for a ± 30% systematic uncertainty, 
we find (see Fig. |^) 

M = 

{asGG/Tt) = 



195J:?5 MeV , 
(330 ± 5 MeV)^ 



-235 ±25 MeV)^ 



in the HV-scheme, and 

M = I95l\l MeV , 
{asGG/n) = (333^^ MeV)'' , 

(qq) = (-245 ± 15 MeV)^ 



(13) 



(14) 



in the NDR-scheme. 

As shown by the light (NDR) and dark (HV) 
curves in Fig. 0, the 75-scheme dependence is con- 
trolled by the value of M, the range of which is 
fixed thereby. The 75 scheme dependence of both 
amplitudes is minimized for Al ~ 190 — 200 MeV. 
The good 75-schcme stability is also enjoyed by 
e'/e and Bk- 

The fit of the amplitude Aq and A2 is obtained 
for values of the quark and gluon condensates 
which are in agreement with those found in other 
approaches, i.e. QCD sum rules and lattice, al- 
though the relation between the gluon conden- 
sate of QCD sum rules and lattice and that of the 
xQM is far from obvious. The value of the con- 
stituent quark mass M is in good agreement with 
that found by fitting radiative kaon decays |14| . 

In Fig. y we present the anatomy of the relevant 
operator contributions to the GP conserving am- 
plitudes. It is worth noticing that, because of the 
NLO enhancement of the 1 = matrix elements 
(mainly due to the chiral loops), the gluon pen- 
guin contribution to Aq amounts to about 20% of 
the amplitude. 
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Figure 2. Test of the 75-scheme stability of the 
Aq and A2 amplitudes as functions of M. The 
light (dark) curves correspond to the NDR (HV) 
results, while the horizontal lines mark the exper- 
imental values. The two figures on the right plot 
AAi = 2|(Af ^ - Af ^«)/(Af^ + Af ^«)|. For 
the values of (qq) and (agGG/Tr) in eqs. (p^)- 
(nj) we find 75-scheme independence for M ~ 
190 - 200 McV. 

AO Ql Q2 03 Q4 Q5 Q6 All 




Ql 02 Q3 Q4 Q5 06 All 

A2 01 Q2 03 Q4 05 Q6 All 

4f 




03 Q4 05 Q6 All 



Figure 3. Operator- by-operator contributions to 
Ao X 10^ and A2 x 10^ (GeV) in the HV scheme. 
In (light) dark the (LO) NLO value. Notice the 
O(p^) enhancement of the gluon penguin operator 
in the / = amplitude. 



Turning now to the A^* = 2 Lagrangian. by 
using the input values found by fitting the A/ = 
1/2 rule wc obtain in both 75-schemcs 



Bk = 1.1 ±0.2. 



(15) 



The result ( [l5[) includes chiral corrections up to 
0{p'^) and it agrees with what we found in |10|| . 

Notice that no estimate of e'/e can be consid- 
ered complete unless it also gives a value for Bk- 
The case of the xQ^i for instance, is telling in- 
sofar as the enhancement of Bq is partially com- 
pensated for by a large B^ (and accordingly a 
smaller Im At ) . 



3.1. The bag parameters Bi 

The bag parameters Bi, defined as 



S, = (0, 



\ model 



/(Q, 



\VSA 



(16) 



have become a standard way of displaying the val- 
ues of the hadronic matrix elements in order to 
compare them among various approaches. How- 
ever they must be used with care because of their 
dependence on the rcnormalization scheme and 
scale, as well as on the choice of the VSA param- 
eters. 

They are given in the xQM in table Q in the 
HV and NDR schemes, at fi = 0.8 GeV, for the 
central value of Aq^j-,. The uncertainty in the 
matrix elements of the penguin operators Q5-S 
arises from the variation of (qq). This affects 
mostly the B^q parameters because of the lead- 
ing linear dependence on {qq) of the Qs^g matrix 
elements in the xQM^i contrasted to the quadratic 
dependence of the corresponding VSA matrix el- 
ements. Accordingly, B^q scale as (qq)^^, or via 
PC AC as mq, and therefore are sensitive to the 
value chosen for these parameters. For this rea- 
son, we have reported the corresponding values 
of ^5,6 when the quark condensate in the VSA 
is fixed to its PCAC value. It should however be 
stressed that such a dependence is not physical 
and is introduced by the arbitrary normalization 
on the VSA result. The estimate of e' is therefore 
almost independent of niq , which only enters the 
NLO corrections and the determination of Bk ■ 

The enhancement of the Qs^e matrix elements 
with respect to the VSA values — the conven- 
tional normalization of the VSA matrix ele- 
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9.3 




9.7 




Bf' 


2.8 




2.9 




Bf ) = Bf' 


0.42 




0.39 




Bi 


-2.3 




-3.0 




Bi 


1.9 




1.3 




Bs -Be 


1.8 X - 
2.6 


-240 McV)3 


1.3 X i 
2.4 


-240 MoV)^ 
(?9> 


o(0) 
-°9 


3.5 




3.4 




n(0) 


4.3 




5.2 




B^'^ ^ Bf ) 


0.89 




0.84 




n(2) _ o(2) 
^9 — ^10 


0.42 




0.39 





Table 1 

Central values of the Bi factors in the HV and 
NDR rcnormalization schemes. For B5 g the lead- 
ing scaling dependence on (qq) is explicitly shown 
for a conventional value of the condensate. All 
other Bi factors are either independent or very 
weakly dependent on (qq). The dependence on 
{agGG/ir) in the ranges of eqs. (p^)-(p^) remains 
always below 10%. 



ments corresponds to taking {qq){0.8 GcV) ~ 
(-220 MeV)3— is mainly due to the NLO chiral 
loop contributions. Such an enhancement, due to 
final state interactions, has been found in 1/A^c 
analyses beyond LO ||l^,|lq], as well as old and 
recent dispersive studies |l7| . 

4. Bounds on ImAt 

The updated measurements for the CKM ele- 
ments iKib/V'cbl implies a change in the determi- 
nation of the Wolfenstein parameter r] that enter 
in Im At . This is of particular relevance because 
it affects proportionally the value of e'/e. 

The allowed values for ImAt — f?|Kis||VcfeP are 
found by imposing the experimental constraints 
for e, iKib/Vcbl, Arrid and Am,,. By using the 
method of Parodi, Roudeau and Stocchi |l|], who 
have run their program starting from our inputs 
listed in fE2|, it is found that 



ImAt = (1.14 ±0.11) X 10""^, 



(17) 



where the error is determined by the Gaussian 
distribution in Fig. H. 
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Figure 4. Bounds on the Wolfcnstcin parameters 
fj = {1-AV2) and p = {l-X'^/2) and distribution 
of Im At according to Parodi et al. for the input 
parameters in p3]. 



5. Estimating e'/e 

The value of e' computed by taking all input 
parameters at their central values is shown in 
Fig. |[ The figure shows the contribution to e' 
of the various operators in two 75 renormaliza- 
tion schemes at /x = 0.8 GeV and 1.0 GeV. The 
advantage of such a histogram is that, contrary 
to the Bj, the size of the individual contributions 
does not depend on some conventional normaliza- 
tion. 

As the histogram makes it clear, the gluon- 
penguin operator Qq dominates and there is very 
little cancellation with the electroweak penguin 
operator Qs- The dominance of the / = com- 
ponents in the xQ^^ originates from the 0{p'^) 
chiral corrections, the detailed size of which is de- 
termined by the fit of the A/ = 1/2 rule. It is a 
nice feature of the approach that the renormaliza- 
tion scheme stability imposed on the CP conserv- 
ing amplitudes is numerically preserved in e'/e. 
The comparison of the two figures shows also the 
remarkable renormalization scale stability of the 
central value once the perturbative running of the 
quark masses and the quark condensate is taken 
into account. 

In what follows, the model-dependent parame- 
ters M, {asCG/n) and (qq) are uniformly varied 
in their given ranges (flat scanning), while the 
others are sampled according to their normal dis- 
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Figure 5. Contribution to e'/e (in units of lO"'^) 
of each penguin operator in the HV and NDR 
75-schemes. The figure on the left (right) corre- 
sponds to /i = 0.8 (1.0) GeV. 



tributions. Values of e'/e found in the HV and 
NDR schemes are included with equal weight. 

For a given set, a distribution is obtained by 
collecting the values of e'/e in bins of a given 
range. This is shown in Fig. n for a particular 
choice of bins. Because the skewness of the dis- 
tribution is less than one, the mean and the stan- 
dard deviation are a good estimate of the central 
value and the dispersion of values around it. 

This statistical analysis yields 



e'/e= (2.2 ±0.8) x 10" 



(18) 



A more conservative estimate of the uncertain- 
ties is obtained via the flat scanning of the input 
parameters, which gives 



0.9 X 10"^ < Re e'/e < 4.8 x 10 



-3 



(19) 



In both estimates a theoretical systematic error of 
±30% is included in the fit of the CP conserving 
amplitudes ^0 and A2. 
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Figure 6. Distribution of values of e'/e (in units 
of 10^"^). Normalized bins are plotted against the 
values of e'/e of each bin. 



Figure 7. Theory vs. experiment in the year 2000. 
The gray band is the average experimental result. 
See m, |0l, ll|, H, lil, H, H, H for 



details on the various estimates. 



The stability of the numerical outcomes is only 
marginally affected by shifts in the value of ^ri+ri' ■ 
Any effective variation of flri+ri' is anti-correlated 
to the value of (agGG/Tr) obtained in the fit of 
A2. We have verified that this affects the calcula- 
tion of Bji and the consequent determination of 
Im At in such a way to compensate numerically in 
e'/e the change of ri,,+,,' . Waiting for a confident 
assessment of the NLO isospin violating effects in 
the K — > TTTT amplitudes, we have used for Qn+ri' 
the 'LO' value 0.25 ±0.10. 

The weak dependence on some poorly-known 
parameters is a welcome outcome of the correla- 
tion among hadronic matrix elements enforced in 
our semi-phenomenological approach by the fit of 
the A/ = 1/2 rule. 

While the xQ^ approach to the hadronic ma- 
trix elements relevant in the computation of e'/e 
has many advantages over other techniques and 
has proved its value in the prediction of what has 
been then found in the experiments, it has a se- 
vere short-coming insofar as the matching scale 
has to be kept low, around 1 GeV and therefore 
the Wilson coefficients have to be run down at 
scales that are at the limit of the applicability of 
the renormalization-group equations. Moreover, 
the matching itself suffers of ambiguities that 
have not be completely solved. For these reasons 
we have insisted all along that the approach is 
semi-phenomenological and that the above short- 



comings are to be absorbed in the values of the 
input parameters on which the fit to the CP con- 
serving amplitudes is based. Because of its sim- 
plicity, the xQM is clearly not the final word and 
it can now been abandoned — as a ladder used to 
climb a wall after we are on the other side — as 
we work for better estimates, in particular, those 
from the lattice simulations. 



6. Other estimates 

Figure M summarizes the present status of the- 
ory versus experiment. In addition to the current 
xQM calculation (and an independent estimate 
similarly based on the xQM), there are eight es- 
timates of e'/e published in the last year. Tri- 
este's, Miinchen's and Roma's ranges are updates 
of their respective older estimates, while the oth- 
ers are altogether new. 

The estimates reported are discussed in P2| . 
The effects first pointed out by the xQM have 
been taken up and refined by some of the other 
approaches, in particular: non-factorizablc cor- 
rections [^Ijlll^,!! and FSI |l|,|l|,||. Fig- 
ure [7| makes clear that, except for estimates heav- 
ily based on VSA for the crucial operator Qe, all 
other estimates are in agreement with the exper- 
iments. 
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